Cooling a single atom in an optical tweezer to its quantum ground state 
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We report cooling of a single neutral atom to its three-dimensional vibrational ground state in an optical 
tweezer. After employing Raman sideband cooling for tens of milliseconds, we measure via sideband spec- 
troscopy a three-dimensional ground-state occupation of ~ 90%. We further observe coherent control of the 
spin and motional state of the trapped atom. Our demonstration shows that an optical tweezer, formed simply 
by a tightly focused beam of light, creates sufficient confinement for efficient sideband cooling. This source of 
ground-state neutral atoms will be instrumental in numerous quantum simulation and logic applications that re- 
quire a versatile platform for storing and manipulating ultracold single neutral atoms. For example, these results 
will improve current optical tweezer experiments studying atom-photon coupling and Rydberg quantum logic 
gates, and could provide new opportunities such as rapid production of single dipolar molecules or quantum 
simulation in tweezer arrays. 



Trapped ultracold neutral atoms provide a promising start- 
ing point for quantum simulation and computation. Ideally, 
experiments would be able to initialize a homogeneous array 
of atoms in an arbitrary quantum state, reconfigure this array 
in real time, and turn on and off interactions. Optical tweezers 
represent an interesting route towards this vision because they 
have the capability to integrate multi-qubit storage, read-out, 
and transport fl]-@). Recently, an idea for near-deterministic 
loading of optical tweezer traps has been demonstrated J5). 
By combining this idea with the capability to detect and cor- 
rect occupation defects, it is a near-term possibility to build an 
ordered neutral-atom array atom by atom [ 6 1 . However, an un- 
realized prerequisite for this vision of a low-entropy tweezer 
array is the ability to place a specific single atom in its vi- 
brational ground state. Such control of single-particle mo- 
tion has been an integral part of trapped-ion experiments for 
decades. In this article, we show that analogous control can 
be realized with a neutral atom. We optically cool a single 
87 Rb atom in a tweezer to its three-dimensional vibrational 
ground state with 90% probability, which is a two order of 
magnitude improvement upon previous experiments employ- 
ing laser-cooled atoms in optical tweezers. 

Motional control of neutral atoms has a rich history, and 
increasingly interest has turned to the problem of single-atom 
control. To date, optical lattices created by standing waves 
of light have been the main platform for realizing motional 
control of collections of single neutral atoms. Approaches 
have included dramatic demonstrations of the Mott insulator 
transition of an evaporatively cooled gas [§), and explo- 
ration of laser cooling collections of atoms in a lattice |j9T fT7l . 
Most recently, in spin-dependent lattices experimenters have 
harnessed microwave signals to cool and control atomic mo- 
tion |fT31 [171 . Our work in an optical tweezer breaks with 
typical lattice experiments, and instead more closely resem- 
bles the sideband cooling and spectroscopy techniques used 
with atomic ions lfl8l . We hold a single atom in a far- 
detuned tweezer trap and apply near-resonant, pulsed cooling 
and spectroscopy light that couples two ground state hyperfine 
levels (Fig.[TJa)). The complete separation of the trapping and 



cooling beams lfl2l [Pfl [161 allows us to realize rapid cooling 
as well as low trap spontaneous emission rates and hence long 
qubit coherence times. 

Already, optical tweezers have been used to realize Ryd- 
berg quantum-logic gates and a variety of protocols coupling 
atoms and photons. In these experiments the thermal motion 
of the atoms has caused deleterious effects, such as dynamic 
light shifts, mitigated atom-photon coupling, and dephasing 
of high fidelity Rydberg gates (3] @] [19] |20| . Attaining full 
three-dimensional motional control ET1 l22l would not only 
strengthen current tweezer applications, but also expand their 
use to experiments that are currently considered only in the 
context of evaporatively cooled gases. For example, one could 
combine two traps and realize significant wavefunction over- 
lap for Feshbach molecule association, and hence creation and 
control of single dipolar molecules ET1 l23ll . Or with suf- 
ficient tunneling between traps and control of a spin degree 
of freedom, one could study strongly correlated physics pro- 
posed in lattices [24] with a platform amenable to studying 
non-equilibrium dynamics. 

The versatility of our tweezer trapping and cooling plat- 
form also presents new opportunities. In particular, a frontier 
in quantum interfaces is realizing strong-coupling of neutral 
atoms to nanopho tonic circuits 125H271 . but trapping and con- 
trolling atoms in nanoscale, near-surface potentials remains a 
considerable challenge. Ground state cooling of an atom in 
an optical tweezer advances this frontier by both providing a 
mobile, highly localized reservoir of single neutral atoms for 
loading nanoscale potentials, and advancing techniques for in- 
situ cooling within these potentials. 

We laser cool to the ground state by employing a technique 
known as Raman sideband cooling [28 29 1. Raman sideband 
cooling operates by iterating on a two step process that ma- 
nipulates the atomic spin and motion, and resolves motional 
transitions even when the trap frequency is less than the elec- 
tronic transition linewidth. In our scheme, we utilize two hy- 
perfine levels of 87 Rb: \F, m F ) = |2, 2) and |1, 1) (Fig.[l}b)), 
where F is the hyperfine angular momentum quantum num- 
ber and m p its projection. The first step is a stimulated Raman 
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FIG. 1. Experimental setup for tweezer trap, detection, and three- 
dimensional motional control, (a) An optical tweezer created with a 
high NA objective lens traps a single neutral atom, and the atom is 
imaged along the z-axis with the same lens. Orthogonal radial axes, 
indicated by x' and y' ', are addressed by Raman beam 1 (RBI) (<r + - 
polarized) and RB2 (7r-polarized), or RBI and RB3 (7r-polarized). 
RB 1 and RB4 (linearly polarized in y-z plane) address the axial di- 
rection. Note we should be able to cool all three axes with a single 
pair of counterpropagating beams, (b) Level diagram for 87 Rb with 
associated beams from (a). The Raman light is 50 GHz red detuned 
of the excited state manifold. Optical pumping (OP) consists of re- 
pump light on the F = 1 — ► 2' transition along with optical pumping 
light onF = 2 -> 2'. 



transition on |2, 2; n) ■<-» |1, 1; n — 1) that reduces the vibra- 
tional state by one quantum while also realizing a spin flip. 
The second step is the dissipative step: The atom is optically 
pumped back to the initial spin-state via a spontaneous Raman 
process, while the photon carries away entropy. By repeating 
these steps, population accumulates in the |2,2;0) state be- 
cause it is dark to the Raman beams and the optical pumping 
light. 

For Raman cooling to be successful, the optical pumping 
step must preserve the reduced vibrational state. The vibra- 
tional excitation probability depends on the motional quan- 
tum number n and the Lamb-Dicke parameter, rj OP = x$k, 
where xq = (fi/2mw) 1 / 2 is the oscillator length for a particle 
of mass m and k is the optical pumping wave number. Raman 
cooling begins with an atom in a mixed thermal state, with a 
temperature corresponding to an average vibrational quantum 
number fi. In this case, the excitation probability due to a sin- 
gle scattered photon scales with (ij^f) 2 = (2n + \)(rf >p ) 2 , 
and hence the Raman cooling efficiency scales inversely with 



Therefore, ground state cooling (achieving n = Q) requires 

(1) low enough initial temperatures before starting Raman 
cooling, and (2) sufficient confinement, i.e. large trap frequen- 
cies. To realize low initial temperatures in our optical tweezer 
trap, we carefully apply the sub-Doppler cooling technique of 
polarization gradient cooling (PGC) QUI . To realize strong 
three-dimensional confinement we use a tightly focused op- 
tical tweezer trap; the trap is formed using far-off-resonant 
light at 852 nm and a 0.6 numerical aperture (NA) objective. 
With 2.8 mW of power in the central focal spot, we measure 
trap frequencies of {u> z , to x /,u} y '}/2n — {30, 154, 150} kHz. 
The large range of frequencies spanned by the axial (z) and 
radial (x, y) dimensions of our trap allow us to explore the 
challenges to robust Raman cooling that accompany increas- 
ing Lamb-Dicke parameters and initial occupations. 

Another challenge specific to the optical tweezer platform 
is effective magnetic fields (vector light shifts) induced by a 
linearly polarized dipole trap with a non-paraxial focus. This 
additional field could (1) dephase field-sensitive transitions, 

(2) disrupt the quantization axis and hence optical pumping 
fidelity, and (3) increase realizable PGC temperatures. For 
our NA, detuning, and typical intensity, we calculate that the 
effective magnetic field is ±0.13 G over ±50 nm, and points 
along the cross-product of the dipole trap axis (z) and its po- 
larization (y) On . During Raman cooling, we set our quan- 
tization field to 3 G in a direction orthogonal to the effective 
magnetic field, mitigating effects (1) and (2). 

Our experiment begins by loading the tweezer trap from a 
magneto-optical trap overlapped with the optical trap focus. 
We use light-assisted collisions to realize zero or one atoms 
in the trap with approximately equal probability JTJ . In order 
to post-select on the presence of one atom, we take an initial 
image of the atom immediately after the loading sequence. 
We start the cooling sequence by applying PGC light for 5 ms 
in a (7 + -(7~ configuration using three retro-reflected beams 
(see the Appendix). We address the challenging problem of 
stabilizing the phase of the interference pattern with respect 
to the sub-wavelength extent of the atom in the tweezer trap 
by modulating the position of our retro-reflecting mirrors at 
1 kHz. This yields a time-averaged cooling that does not ex- 
hibit shot-to-shot variations in the final PGC temperature due 
to slow fluctuations in the interference pattern. 

As our first application of Raman coupling, we diagnose a 
thermal occupation using sideband spectroscopy. After PGC 
is performed, we optically pump to \F,mp) = |2, 2). A sin- 
gle pair of Raman beams is chosen according to the dimension 
we wish to probe (Fig.[TJa)), and we interrogate the atom with 
a square pulse of length At, leading to transitions to the |1, 1) 
state when a resonance is satisfied. The Lamb-Dicke param- 
eter for the Raman transition, which quantifies the motional 
coupling, is given by r) R = x^Ak, where Afc is the momen- 
tum transferred by the Raman process along the motional axis. 
For our geometry rj R = 0.22 for the radial dimension and 
i] R = 0.23 for the axial dimension. 

Figure|2|a) shows sideband spectra after PGC along orthog- 
onal radial directions. The asymmetry between the An = — 1 
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FIG. 2. Single atom sideband spectra and Rabi oscillations in the radial dimensions before (a, b) and after (c, d) ground state cooling, (a) 
The black squares are a carrier peak in the y direction using a At = 15 us (near n) pulse. The red circles (orange triangles) are sidebands 
along the y (x 1 ) axis for a 75 fis (near it) pulse, demonstrating an initial thermal population of vibrational states. The solid lines are fits to a 
Rabi sine function (See Supplementary Material for definition) for the carrier and Lorentzians (an approximation) for the sidebands; each fit 
contains an offset at our measured background (gray shaded region centered at 0.04). (b) Carrier Rabi oscillations for the y' direction showing 
dephasing of a thermal state. Here the carrier Rabi frequency was set to 15 kHz, instead of the usual 26 kHz. The solid line is a fit to the data 
using a thermal distribution of Rabi frequencies, (c) Raman cooled radial sidebands; no Raman cooling is applied to the axial direction for 
these data. The black squares are a cooled carrier peak using a 15 /j,s pulse. The blue circles (green triangles) are spectra along the y' (x') 
axis using a 75 /is pulse, displaying a significant asymmetry that is the hallmark of a large ground state population, (d) Rabi oscillations for 
a radial ground state cooled atom with a fit to a damped sine for the carrier (black squares) and the An = +1 sideband (blue circles), which 
demonstrates coherent control of the spin-motional states; the carrier dephasing is suppressed due to the purity of the vibrational distribution. 
Each data point is an average of 150 experimental runs, and hence ~ 75 atoms. 



and the An = +1 peaks determines the population of the 
atoms in the ground state. (Because we begin in the upper hy- 
perfine state, the An = —1 peak is on the right in Fig. [2]) By 
equating the ratio of the An — — 1 and An = +1 sidebands 
to n/(n + 1) we find n y > = 1.1 ± 0.4 or T> = 11 ± 3 uK 
assuming a thermal population distribution; correspondingly 
n x i = 1.7 ± 0.7. Using the dephasing of the carrier Rabi 
oscillations in these data, we can also extract a temperature 
via the coherent evolution of the thermal state [Fig. |2|b)] to 
find T y / < 16 ± 2 uK l32l . For comparison, we also em- 
ploy a standard thermometry technique in which the atom is 
quickly released, and the probability of recapturing the atom 
at a variable time later is measured and compared to a classical 
Monte Carlo model [33 1 from which we estimate a tempera- 
ture of 13 ± 1 uK. The agreement between our three mea- 
surements validates sideband spectroscopy as a reliable form 
of thermometry in an optical tweezer, and we find low PGC 
temperatures are possible despite the varying effective mag- 
netic fields within the focus. 

For the radial dimensions, the PGC allows us to start our 
Raman cooling reasonably far into the Lamb-Dicke regime 
with r/Off ~ 0.3 (n ~ 1.5 and n OP = 0.16). Like ion ex- 
periments, our cooling scheme operates in a pulsed format: 



We apply a coherent Raman transition on the An = — 1 side- 
band using an intensity corresponding to a ground state carrier 
Rabi frequency of £l c = 2ir x 26 kHz, where the An = — 1 
sideband Rabi frequency is il s i, ~ n R n c y/n. Therefore, we 
apply 60 us pulses for 47 cycles and switch to 75 {is pulses 
for the last 3 cycles when the remaining excited state fraction 
is primarily in the first excited state. The pulses are applied 
to each radial axis by alternating Raman beam pairs. Inter- 
spersed between the Raman pulses are optical pumping pulses 
90 us long that recycle the atom back to |2, 2). Note we do 
observe cooling of both radial dimensions even if we employ 
a single pair of cooling beams, indicating there is coupling 
between the radial dimensions of our anharmonic potential. 
To assure and verify cooling in both dimensions, however, we 
cool and probe each axis separately. 

After Raman cooling we see a significant asymmetry in 
the radial sideband spectroscopy due to a large ground state 
occupation (Fig. |2jc)). While the An = —1 sideband is 
suppressed, the An = +1 sideband has increased in height 
due to decreased dephasing as the thermal distribution is nar- 
rowed. Figure |2|d) shows Rabi oscillations for both the car- 
rier (black) and An = +1 sideband (blue) transition, the 
latter of which oscillates slower by a factor of n R . Fig- 
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FIG. 3. Axial spectra and 3D ground state cooling of a single neutral atom, (a) A thermal axial spectrum (red squares) using an intensity 
corresponding to a cold carrier Rabi frequency of 12 kHz and a At — 65 fis Raman pulse. The data are fit to Lorentzians (solid line) to guide 
the eye. (b) Result of simultaneous sideband cooling in three-dimensions, demonstrating significant sideband asymmetries and simultaneous 
ground state occupations in all dimensions. The axial data (center) illustrates our cooling parameters, and is performed with a carrier Rabi 
frequency of 10.6 kHz and a pulse of 236 /is (near 7r pulse on the ground state An — +1 sideband, a 5n pulse on the carrier). For the radial 
data, the blue circles (green triangles) are spectra along the y' (x 1 ) axis using a 75 pis pulse. The solid lines on the An = +1 sidebands 
are Lorentzian fits, (c) After 3D cooling, axial spectroscopy for a halved carrier Rabi frequency of 5 kHz and a pulse of 450 as. Better 
spectroscopic resolution affirms a large axial ground state occupation. 



ure |2jd) further demonstrates the coherence of our motional 
transitions. The carrier decays less quickly than in Fig. |2jb), 
as expected for colder atoms; the sideband transition decays 
slightly faster than the cold carrier due to its relative nar- 
rowness 1341 , but note that it maintains high contrast on the 
first oscillation. To assess our final occupation, we compare 
the measured signal at the position of the An = — 1 transi- 
tion peak to the measured background level due to any atoms 
left in F = 1 and imperfect push-out efficiency. We find 
{^,fv} = {0.05tg;S| ) 0.02l«;g*}. 

To achieve large three-dimensional ground state occupa- 
tions, we must cool the weaker axial dimension of the trap 
where both spectroscopy and cooling are more challenging 
due to the smaller trap frequency. Figure |3ja) shows a ther- 
mal axial mode spectrum after PGC. The near equality of the 
An = +1 and An = — 1 transitions and the presence of sig- 
nificant second order sidebands (given rj R = 0.23 for this di- 
mension) are consistent with a small initial ground state pop- 
ulation. Assuming an isotropic initial temperature of 12 uK, 
we would expect n — 8, which corresponds to a challenging 
starting point of n^ p = 1.4. Further, the smaller trap fre- 
quency makes it difficult to spectroscopically separate the car- 
rier and sideband peaks while maintaining Rabi frequencies 
that are insensitive to technical dephasing. For the cooling, a 
large Rabi frequency leads to off-resonant carrier transitions 
that cause heating, while too small a Rabi frequency leads to 
smaller transfer efficiencies and slow cooling. 

Despite these barriers, we are able to Raman cool in the 
axial dimension and achieve the significant three-dimensional 
ground state occupations evidenced in Fig. [3jb) after 33 ms 
of cooling (see the Appendix for pulse parameters). Unlike 
the radial directions, we use different parameters for the ax- 
ial cooling than are ideal for spectroscopy. To highlight this 
distinction, we show two spectra of cooled atoms. The axial 
mode spectrum in Fig. |3jb) approximates the parameters used 
during our cooling, which balance speed with spectroscopic 



resolution. While there is a clear sideband asymmetry in the 
spectrum, the size of the Rabi frequency compared to the ax- 
ial trap frequency complicates a temperature analysis because 
off -resonant carrier transitions occur at the frequency position 
of the sidebands. This is illustrated in Ref. [35] where we 
show calculated spectra for these pulse parameters, but a full 
comparison to this calculation is complicated by the partial 
dephasing we observe. 

In Fig. [3j c X we halve the spectroscopy Rabi frequency to 
sacrifice coherence for spectroscopic resolution. For this Rabi 
frequency and pulse, the spectrum is dephased and hence 
we can understand the spectrum simply as a set of multi- 
ple Lorentzians. We extract our axial temperature from the 
sideband asymmetry in Fig[3jc), assuming the dephasing uni- 
formly affects the two sideband peaks. For both the radial and 
axial dimensions we analyze the data attributing all of the sig- 
nal observed at the An = — 1 position to the sideband, thus 
placing an upper bound on the achieved temperature. Tak- 
ing this approach we extract temperatures in all three dimen- 
sions, again subtracting off a measured background and as- 
suming a thermal distribution. The result is {n x > , n y > , n z } = 
rn no+0.07 q 01 +o.06 q 08+ - 08 ; 

\U.UZ_ Q2 , U.U-L_q 01 ,U.UO_ QgJ-. 

The above occupations indicate we have cooled a single 
neutral atom to the ground state of an optical tweezer with 
97^% probability in the radial plane, and 93^% proba- 
bility in the axial direction, and hence a three-dimensional 
ground state population of 90j^ 6 %. While this estimate ne- 
glects imperfections in our spin preparation and detection, this 
is not a fundamental impediment to the final temperature. The 
current limitation to our quoted temperature is the precision 
of the spectroscopy measurement. Investigation of tunneling 
between adjacent tweezer traps would allow qualitatively new 
ways of analyzing the purity of the ground state preparation, 
because bosonic enhancement of tunneling rates is contingent 
upon this preparation. Importantly, we have measured an up- 
per bound to the heating rate of less than 1 radial vibrational 
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quanta per second, suggesting tunneling experiments are fea- 
sible with sufficiently tight, adjacent tweezer traps. Such ex- 
periments would open a new avenue of research for ultracold 
atoms in optical tweezers, and establish their utility for several 
directions in quantum information and bottom-up approaches 
to quantum simulation. 
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Note added - Related concurrent studies independently car- 
ried out by another group are described in Ref. l37l 

APPENDIX 

Optical tweezer trap and single atom loading. The optical 
tweezer is formed by a custom high-numerical objective lens 
that is designed to be diffraction limited at both 850 nm and 
780 nm ll36l . From our measured trap frequencies, we infer 
a 1/e 2 beam radius of ~ 0.76 /im and a depth of 1.4 mK for 
our intensity. Abberations dominated by astigmatism slightly 
increase our spot size compared to a diffraction limited spot 
for a 0.6 NA lens. We measure a 1/e lifetime of a single atom 
of about 5 s, which can be extended at the expense of our trap 
load rate. 

We load the tweezer trap from a vapor cell magneto-optical 
trap (MOT) that has filled for 100 - 200 ms. The MOT con- 
sists of two orthogonal pairs of retro-reflected beams along x' 
and y', as well as a third pair at approximately 45° to the z- 
axis and 10° above the y — z plane. For polarization gradient 
cooling (PGC), the same beams are used as for the MOT light 
and the bias fields are set for zero magnetic field; we zero our 
fields on the basis of microwave spectroscopy. The PGC de- 
tunings are optimized for each stage of the experiment. For 
loading into the tweezer trap, the F = 2 — 3' light is detuned 
60 MHz red of the bare optical transition; the PGC to prepare 
for Raman cooling uses light 23 MHz red of the bare optical 
transition [ 33 1 . 

All of our data rely on the ability to make spin-sensitive 
measurements of our 87 Rb atom. Our procedure is to apply 
light resonant with the light-shifted F = 2 — > 3' transition 
while the atom is in the trap to retain atoms only in F = 1. 
We then image the F = 1 atoms by applying PGC (including 
repump) light. This allows us to collect fluorescence for 25 to 
50 ms while imaging this light onto an CCD camera with ~ 
7% overall efficiency with the same high-NA lens that creates 
the tweezer trap. 

3D cooling parameters. For the three-dimensional cooling 
shown in Fig. |3jb),(c) we toggle the Raman laser pulses to 
address each of the axes of the trap, and between each pulse 
we insert 90 [is of optical pumping. The cooling process in 
total occurs in 75 cycles: The first 50 cycles use a Raman 



pulselength of At — 48 [is (At = 40 /is) for the radial (axial) 
directions, and then At = 72.5 /is (At = 80 /xs) for the final 
25 cycles. We use an intensity corresponding to a cold carrier 
Rabi frequency of 13 kHz for the axial direction, and 31 kHz 
for the radial dimension. 
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SUPPLEMENTARY MATERIAL 

Optical layout. Figure SI illustrates the basic optical 
layout of our experiment. The optical tweezer that pro- 
vides three-dimensional confinement is created by focus- 
ing 852 nm light through a custom built objective lens. 
It has a working distance of 21 mm, and is corrected for 
the 0.25 in. fused-silica window of our octagonal vacuum 
cell. The objective consists of both aspheric and spheri- 
cal elements. In creating the light for the trap, we incor- 
porate two acoustic-optic modulators (AOMs) that yield 
deflections in two orthogonal directions. This will allow 
for future dynamic positioning of multiple traps and in- 
dividual depth control of these traps in real time. The 
same objective that creates the traps also collects fluo- 
rescence from the trapped atom, which is then imaged 
on to a CCD array. 
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FIG. SI. Optical layout for the optical tweezer trap and atom 
detection. 



Raman transitions. The Raman beams are detuned 
approximately 50 GHz red of the D2 line with an ex- 
ternal cavity diode laser, which is passively stable. The 
intensity of each of the four Raman beams is separately 
switched with acoustic optic modulators, and RBI ad- 
ditionally passes through an in-fibcr electro-optic mod- 
ulator (EOM) operating near the hyperfine splitting of 
6.8 GHz. RBI, prior to passing through the EOM, is fre- 
quency shifted by +10 MHz from RB2, RB3, and RB4. 
This ensures that only one of the optical sidebands on 
RBI is resonant with a Raman process while the other 
optical sideband is 20 MHz off-resonance, thereby pre- 
venting quantum interference effects and transitions as a 
result of any residual <j~ polarization on RBI [1]. The 
10 MHz frequency difference between the two carriers is 
also sufficiently large so as to avoid parametric heating. 
As discussed in the text, RBI is always used in combina- 
tion with one of Raman beams to address an axis. The 
beam detunings and polarizations are such that start- 
ing in F — 2 the atom always stimulated emits into 
RBI while absorbing a photon from the other beam in 
use, which determines the direction of momentum trans- 
fer effected with each beam pair. The light shift on the 
|1, 1) O |2, 2) transition is ~ 90 kHz for all of the Raman 
beam pairs, which is dominated by the vector light shift 
induced by the cr + -polarized RBI beam. 

When driving Raman transitions we observe some 
technical dephasing; we find the dephasing rates increase 
as the Rabi frequency of the driven transition is reduced. 
This is seen in the different dephasing rates of the car- 
rier and the An = +1 Rabi oscillations in Fig. 2(d), as 
well as the dephasing observed in the axial spectra. We 
suspect decreasing the resonance width increases sensi- 
tivity to technical variations in the light-shifted hyper- 
fine resonance. This instability could be attributed to 
intensity fluctuations in the Raman beams or the optical 
tweezer vector light shift, which varies across the atomic 
wavepacket. 



Optical pumping setup and optimization. The 

heating during the optical pumping step determines the 
efficiency of Raman cooling. In our case, this heating 
rate depends on the recoil per photon and the number 
of photons scattered to optically pump the atom from 
\F,rriF) = |lj 1) to |2, 2). Polarization impurity increases 
the number of photons scattered, and diminishes the fi- 
delity of the pumping. 

We take the following steps to optimize the polariza- 
tion seen by the atoms. The polarization is set with a 
polarizing beam splitter and a quarter-wave plate. Af- 
ter aligning the optical pumping beam to the atoms, we 
optimize the direction of the quantization axis to be par- 
allel to the direction of the optical pumping beam. We 
first optically pump to |2, 2) using F = 2 — 2' light and 
repumping F = 1 — 2' light. Subsequently, we apply the 
same F = 2 — 2' light in the absence of the repumping 
F = 1 - 2' light. If |2,2) were truly dark, this would 
have no effect. However, we can detect residual tt and 
a~ light via the depumping rate out of |2,2). Since the 
F = 2' level decays with equal probability to F = 1 and 
F = 2, we can detect depumping via population build 
up in F = 1. As shown in Fig. S2, we then optimize 
the direction (field angle) and magnitude (tip field) of a 
small transverse field added to the 3 G quantization axis 
along the y-axis by minimizing this depumping rate. 

Calculation of optical pumping photons scattered. 

To calculate the feasibility of the number of cooling cy- 
cles required in the experiment, we determine theoreti- 
cally the number of optical pumping and repumping pho- 
tons scattered in the ideal situation of pure, resonant a + 
F = 2-2' light and repumping F =1-2' light [2]. For 
this calculation, the Hilbcrt space consists of the mag- 
netic sublevels of F = 1, F = 2, and F = 2'. We per- 
form a master equation calculation, given the atomic line 
strengths, by numerically solving the differential equa- 
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FIG. S2. Optical pumping optimization. After optical pump- 
ing to 1 2, 2), we vary the magnetic field quantization axis 
transverse angle (a) and magnitude (b) and observe the 
depumped fraction after 4 ms of exposure to F — 2 — 2' light. 
We attribute residual depumping to imperfect quarter-wave 
plate alignment, and finite detuning from other transitions on 
the D2 line. 



tion defined by a Lindblad operator C, 
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where p is a density operator describing the sub-level 
populations, the 7.^ are defined as, 
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and Tij describes the decay rate from \i) to \j) due to 
spontaneous emission, or the incoherent scattering rate 
of \i) to \j) via the optical pumping and repump beams. 
From the time evolved density matrix, the number of 
scattered photons (n 7 ) can be determined from 
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where the sum i is taken over the excited state manifold 
F = 2' and the sum j is taken over the sublevels of the 
two hyperfine ground states F = 1 and F — 2. A plot of 
photon number scattered and population in |2, 2) versus 
time is shown in Fig. S3. We emphasize that (n 7 ) is 
minimized when the repump light is also cr + -polarized. 

Spectroscopy lineshape For the carrier spectra shown 
in Fig. 2(a) and Fig. 2(c) we fit to the expected the- 
oretical lineshape. For a two level system driven by a 
coherent source with Rabi frequency Q, pulse At, and 
detuning 5, the transition lineshape is given by the Rabi 
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FIG. S3. Optical pumping calculation. Numerical results for 
optical pumping rate into |2, 2) and the number of photons 
scattered as a function of time. The parameters for this cal- 
culation are a repumping rate of 1 MHz, optical pumping rate 
of 100 kHz, and excited state linewidth of 6 MHz. 



sine function, 
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where A describes an overall amplitude. For the thermal 
carrier spectra displayed in Fig. 2(a) this model is an 
approximation because it assumes a pure Rabi frequency. 
In Fig. 2(c) the Rabi frequency is primarily that of the 
ground state carrier transition and hence this model is a 
good description. 

Simulation of coherent dynamics under Raman 
coupling interaction. We determine the coherent dy- 
namics in order to better understand our measured spec- 
tra and cooling efficiency. We approximate the internal 
degrees of freedom as a spin- 1 /^ system, and the exter- 
nal degree of freedom we truncate at 10 vibrational lev- 
els. For our initial temperatures in the radial degrees of 
freedom, the highest level has a population of < 0.005, 
and for the axial < 0.05. Our calculation of spectra and 
dynamics only considers one external degree of freedom 
(ID), where the dimension of interest is distinguished by 
its associated trap frequency. In the rotating frame, the 
Raman interaction is written, 
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where f2 is the bare Rabi frequency, a + is the spin-!/2 
raising operator, AA; is the momentum transfer for the 
Raman beam configuration, and the exponential term is 
the momentum translation operator associated with the 
coherent momentum kick due to absorption of a photon 
from one beam and stimulated emission into the other [3] . 
The argument of the exponential can be recast as follows, 



iAki = iry fi (a t + a), 
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the ground state wave-packet size, and a (a') the 
annihilation (creation) operators of vibrational excita- 
tions. For numerical calculations, we Taylor expand the 
momentum-translation operator in powers of r) R to 4th 
order. In the rotating frame, the bare Hamiltonian is, 
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We solve numerically the combined Hamiltonian H 
Hq + iJpt, such that for a density matrix p, 



p(t) = -i[H,p], 



(S10) 



for variable initial conditions. For example, we can solve 
for the evolution of thermal states, 
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where Z is the partition function associated with a trun- 
cation of 10 states. 

We can numerically simulate spectra and dynamics 
with this model. We show in Fig. S4(a) the expected 
spectrum for an axial thermal distribution characterized 
by n = 0.1 and the spectroscopy parameters used in the 
spectrum in Fig. 3(b); this calculation does not include 
effects of technical dephasing. For these pulse parame- 
ters, we see significant off-resonant carrier contributions 
for our square pulses, which are comparable to the height 
of the An = — 1 sideband. This behavior supports our 
discussion in the text, namely that carrier contributions 
can complicate a temperature analysis. 

Now we consider the analysis of the thermal dephas- 
ing observed in carrier Rabi oscillations as in Fig. 2(b). 
While it is possible to use Eqn. S10 to numerically calcu- 
late carrier and Rabi evolutions, there are also exact ex- 
pressions for the Rabi frequencies. We use these expres- 
sions to extract a temperature from carrier Rabi evolu- 
tions. For a thermal distribution of populated vibrational 
levels, there is dephasing of the carrier Rabi oscillations 
because each level has a different associated carrier Rabi 
frequency to above first order in rj R . To fit thermal car- 
rier Rabi oscillation data such as that shown in the main 
text in Fig. 2(b), we use the following model with fio, T, 
and an overall normalization as free parameters, 
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FIG. S4. Coherent dynamics, (a) We simulate an axial 
spectrum for the parameters: n = 0.1, Qo/2n = 10.6 kHz, 
At = 236 ps, r) R = 0.23, and lo/2tt = 30 kHz. (b) For the ra- 
dial direction, which exhibits r] R — 0.22 and u)/2ty = 150 kHz, 
we show the carrier Rabi evolution for temperatures of 1 pK, 
10 pK, and 20 pK corresponding to the blue, green, and red 
plots, respectively. 



where the level dependent carrier Rabi frequency can be 
expressed, 



fi(n) = r2 oe -^) 2 / 2 L„(( ?7 «) 2 ), 



(S13) 



where L n is the Laguerre polynomial, and Pt(ji) is the 
normalized Boltzmann factor for a thermal state of tem- 
perature T . Three examples of Rabi evolutions for dif- 
ferent thermal occupations are displayed in Fig. S4(b), 
which demonstrate increasing dephasing as the thermal 
distribution broadens [1, 3-5]. 

Master equation calculation of Raman cooling. To 

understand the expected cooling rate, we model the cool- 
ing process with a master equation formalism. According 
to our optical pumping calculations, we assume there are 
three photons scattered during the optical pumping step 
of the Raman cooling. Hence, we define an optical pump- 
ing operator as follows. We take the cooling state to be 
the spin down state, where Raman transitions cause spin- 
flips to the spin up state. The optical pumping operator 
for an effective optical pumping rate T is defined as fol- 
lows, 



6 = r^V 
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where the spin-lowering operator a recycles the spin- 
state while the second term realizes the associated heat- 
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FIG. S5. Cooling calculation. Starting from a temperature of 
15 /iK, we simulate the cooling process. We assume here that 
three photons are scattered during the optical pumping, and 
that the spontaneous emission process is ignored (only the 
absorption is accounted for). This corresponds to a 150 kHz 
trap frequency, as in our radial direction, and a £2o/(27r) = 
26 kHz Rabi frequency. The coherent pulse lengths on the 
An = — 1 sideband are initially 60 ps, and then 75 pis for the 
last 10 cycles. 



ing due to three scattered photons of wave-number k. We 
accordingly define the Lindblad operator, 

C = OpO^ - \&6p - \p&6- (S15) 
The master equation to be numerically solved is then, 
p{t) = C. (S16) 



To simulate the cooling, we alternate between evolution 
under the coherent process defined by Eqn. S10 and the 
dissipative process denned by Eqn. S16. This accurately 
models our experiment where the optical pumping and 
coherent Raman drive are never concurrent, in contrast 
to a steady state cooling process in which coherent cou- 
pling and optical pumping are simultaneous. We can 
calculate the ground state fraction by finding the expec- 
tation value of the dark state projector. An example of a 
cooling trajectory similar to our experiment is displayed 
in Fig. S5, from which we conclude that the 50 cooling 
cycles we use experimentally is reasonable. With better 
choice of final pulselengths it is possible to increasingly 
improve the dark state population in simulation, and the 
disparity in final occupations between data and exper- 
iment could be attributed to small differences between 
the experimental and calculation parameters. 
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